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■^ 1 INTRODUCTION 



ABSTRACT 

Recently, gamma-ray telescopes AGILE and Fermi observed several middle-aged su- 
pernova remnants (SNRs) interacting with molecular clouds. It is likely that their 
gamma rays arise from the decay of neutral pions produced by the inelastic collision 
between cosmic rays (CRs) and nucleons, which suggests that SNRs make the bulk 
of Galactic CRs. In this paper, we provide the analytical solution of the distribution 
of CRs that have escaped from a finite-size region, which naturally explains observed 
broken power-law spectra of the middle-aged SNRs. In addition, the typical value of 
the break energy of the gamma-ray spectrum, 1-10 GeV, is naturally explained from 
the fact that the stellar wind dynamics shows the separation between the molecular 
clouds and the explosion center of about 10 pc. We find that a runaway-CR spectrum 
of the four middle-aged SNRs (W51C, W28, W44 and IC 443) interacting with molec- 
ular clouds could be the same, even though it leads to different gamma-ray spectra. 
This result is consistent with that of recent studies of the Galactic CR propagation, 
and supports that SNRs are indeed the sources of Galactic CRs. 

Key words: acceleration of particles - cosmic rays - ISM: individual objects (W51C, 
W28, W44, IC 443) - shock waves - supernova remnants. 
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Supernova remnants (SNRs) 
Galactic cosmic rays (CRs). 



are thought as the origin of 
The most popular accelera- 



tion mechan ism at SNR is the diffusive shock accelera- 
tion (PSA) (lAxford et al.lll977l : lKrvmskvl[l977l : lBeilll978l : 
iBlandford fc OstrikeiJIigTa ). To make the Galactic CRs, a 
few tens of percent of the ex plosion energy of a sup ernova is 
converted to the CR energy (jBaade &: Zwickvll 19341 ). Recent 
gamma-ray telescopes, such as Fermi and AGILE, suggest 
that middle-aged SNRs interacting with molecular clouds 
likely emit hadronic gamma-rays, which arise from the decay 
of tt" produced in inelastic c ollisions of the accelerated pro- 
tons with molecular clo uds ("Abdo et aJ] l2009bl . I2010a|[bll3 : 
iTavani et aLll201(]| : [GiuU ani et al. 20l3^ This fact supports 
that the SNRs make the bulk of the Galactic CRs. Fermi 
observations revealed that the gamma-ray spectra of the 
middle-aged SNRs have a broken power law form, that the 
break photon energies are typically ~1-10 GeV, and that the 
spectral indices above the break energy are different from 
each other. 
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Several theoretical interpretations of the gamma-ray 
observatio ns of middle-age d SN Rs have been pro posed. 
Aharonia n fc AtovanI (| 19961 ) and iGabici et"all ( |200S^ ) have 
investigated the gamma-ray spectrum from molecular clouds 
illuminated by CRs which have escaped from a nearby SNR. 
They showed that the spectrum is steeper than predicted 
by DSA at the shock bec ause of th e energy dependent 
diffusion. In addition, Aharonian fc At ovan ( 199g) showed 
that the spectrum has a broken power law form due to 
the finite size of th e emission region. On the other hand, 
iMalkov et al.l l|2010l ) recently proposed a possible explana- 
tion. They claimed that in the partially ionized medium, 
CRs accelerated at the shock has a broken power-law spec- 
trum — above the break energy, the damping of waves that 
resonantly scatter the CRs is significant, making the steeper 
spectrum than below the break energy. 

In this paper, we propose another explanation for 
the observed middle-aged SNRs interacting with molecu- 
lar clouds, by reinvestigating the distribution of CRs that 
have escaped from an SNR. So far, the point source ap- 
proximation has been widely adopted. However, for SNRs 
interacting with molecular clouds, the extension of the CR 
sources is important. In section 2, taking into account the 



Y. Ohira et al. 



effect of a finite-size source, we first provide the distribution 
of CRs that have escaped from a SNR. Then, we consider 
the dynamical evolution of CRs interacting with molecular 
clouds, because the CR spectrum that has escaped from an 
SNR depends on t he time evolution of the maximum energy 
ijOhira et al.ll2010l ). In section 3 and section 4, we provide 
two models of the evolution of the maximum energy. The 
gamma-ray spectra from molecular clouds are calculated in 
section 5. Section 6 is devoted to a discussion. 



2 DISTRIBUTION OF CRS ESCAPING FROM 
A SNR 

In this section, we derive the distribution of the runaway 
CRs at a given distance r from the SNR center and at a given 
time t from the explosion time, f{t,r,p), where p is the CR 
momentum. We here assume that the system is spherically 
symmetric and the diffusion coefficient in the interstellar 
medium -Dism(p) is spatially uniform and depends on the 
CR momentum. Then, we solve the diffusion equation given 
by 



9/ 
dt 



{t,r,p) - DisM(p)Af{t,r,p) = qs{t,r,p) 



(1) 



where q ^ is the source term of CRs. Considering the es cape 
process jPtuskin fc Zirakashvilill2005l : IOhira et al.ll2010l '). the 
CRs with a momentum p escape from an SNR at f = tesc(p). 
The Green function, that is, the soluti on in the point sour ce 
case, qs = Ncsc{p)5{r)5{t — tcsc(p)), is (|Atovan et al.lll995l 'l 



/point (i,r,p) = 



where 



^Rd(*.p)y 



7r3/2i?d(t,p)3 



iVcsc(p) 



-Rd = v/4AsM(p)(t-icsc(p)) 
and 



A^osc(p) 



dt d r qs{t,r,p) 



(2) 



(3) 



(4) 



Now, taking into account the fact th at the CRs escape from 
the SNR surface (jOhira et al.ll2010l ). the source term is re- 
placed with 

'7= = ^^^<5(r--i?c.c(p))5(t-tesc(p)) , (5) 

where -Rcsc(p) is the radius where the CRs with momentum 
p escape from the SNR. Then, we find the solution to equa- 
tion IT]) as 



/cxt (t,r,p) 



d^r7pomt(t, |r- r'|,p) 



5(r'-i?esc(p)) 



4jrr' 



Hd(*.P) J 



-RdC.P) y 



47r3/27?d(t,p)i?csc(p)r 



-iVe.c(p)(6.) 



Note that /cxt has a spectral break at p = pbr.cxt because 
/oxt ~ Ncsc/4:n^^'^ RdRcacT for r ~ Rcsc > Rd, otherwise 
/cxt ~ /point, where pbr.oxt is obtained from 



-Rose (pb 



— -Rd(Pbr,cxt) 



(7) 



Below Pbr.cxt (^csc > ^d), the effect of a finite-size source 
becomes important. 



The observed gamma-ray spectrum should be calcu- 
lated by the volume-integrated CR spectrum. In the case 
of the SNR-molecular-cloud interacting system, the emis- 
sion region is a dense cloud or cavity wall, whose density is 
high in the region Li ^ r 5C L2. 

Considering the stellar wind before the supernova ex- 
plosion, the inn er radius of the mole cular cloud Li is about a 
few tens of pc (We aver et al.lll977l ). The volume-integrated 
spectrum of CRs in the region Li ^ r ^ L2 at the SNR age 
fagc is calculated as 



F(p) 



/cxt(tagc,r-,p)47rr dr 



(8) 



Substituting Eq. ((6)) into the above, we obtain 



2 1 y^i?csc 

'\ «d ) 



+ e 



+erf (■ 



Re: 



Rd 
L2 + Re 



O-Ktt^) 



Rd 
Li + Re 



Rd 



0-»'(^4f-)} 



(9) 



where erf (a;) — {2/^/tt) J e~^ dy is the error function. 

There are four characteristic momentum regimes deter- 
mined by the comparison between the distances that CRs 
can reach, iZcsc(p) -f- -Rd(p) and I/i,2(> i?esc(p)), and by the 
comparison between Rd and i?csc; 



(i) 1/1,2 < ReBc{p) + -Rd(p) and 7?esc < Rd, 

(u) Ll,2 < -Reac(p) + -Rd (p) and Resc > Rd (il 

(iii) Li < -Rcsc(p) + -Rd(p) < L2, 
(iv) and Rcsc{p) + Rd{p) < ii,2. 



^Cl 



Let us define Pbr,2 and pcut as the solutions to the following 
equations: 



L2 = i?csc(Pbr,2) + -Rd(Pbr,2) 
Ll = i?csc(Pcut) + -Rd(Pcut) 



(10) 
(11) 



We may expect pcut < Pbr,2 as long as the momentum de- 
pendence of Rcsc is weak enough and/or Rd is more impor- 
tant. Then, a significant fraction of CRs with energies above 
Pbr,2 leave the molecular cloud, while CRs below pcut do not 
essentially reach the cloud. 

Furthermore, we approximate e~^ ~ 1 — x^ + x'^/2 and 



erf(x)Ri(2/0F)(j;- 
erf (a;) « 1 for 1 < a; 


- x'^S) for < 1 < 1, and e""= r; and 
. Then, for Pbr.cxt > Pbr,2, we obtain 


F{p) OC iVcsc(p) X ■ 


^d(p)-^ 
Rd{p)-'B 

P° 



(P > Pbr.cxt) 
!csc(p)"^ (Pbr,cxt > P > Pbr,2 
(Pbr,2 >P >Pcut) 
(Pcut > p) 


For pbr,2 > Pbr.cxt, wc obtaiu 




F{p) OC iVcsc(p) X - 


Rd{p)-' 

p" 




(P > Pbr,2) 

(Pbr,2 > P > Pent) . (13) 

(Pcut > p) 



Note that there is no break at pbr.cxt since the effect of a 
finite-size source is smeared by spatial integration. 

The spectral breaks at Pbr.cxt and Pbr,2 come from the 
finiteness of the source and the emission regions, respec- 



(12) 
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tively. To estimate the values of Pbr,2 and Pcut, we first as- 
sume the diffusion coefficient of the interstellar medium 

Asm(p) = fO^« X f TT^V cm^ «"' . (14) 



Hence we assume that CRs with the knee energy, cpknoc ~ 
]^q15-5 gY^ escape at f = fsodov We adopt a phenomenologi- 
cal approach based on the power-law dependence, 



where x is constant. CR propagation mo dels require Y ~ f 
and S ^ 0.5 as the galactic average value (IBerezinskii et al.l 
[l990|). Assuming that Resc and iesc are constant with p (be- 
cause p-dependence of _Resc is weak), the typical value of pcut 
is estimated as 



Pcsc — Pkn 



/_Rsh_y 

V -Rscdov / 



(18) 



where -Rsodov is the SNR radius at f = iscdov We assume 
a — 6. 5 in order to Pesc = 1 GeV/c at the end of the Sedov 
phase (|Ohira et al.ll20ld '). From equations. (fT7|) and ([T5)). we 
obtain 



Pcut = 7 



(f) 



-l/S 



Ll — Re: 

5 pc 



2/S 



f* f'CSC 

10* yr 



-l/S 



GeV/c , (15) Do = 



where we assume 5 — 0.5, and Pbr,2 is obtained by replacing 
Ll with L2. Although Li is order of 10 pc, AGILE observa- 
tions show that a fewx 100 MeV photons come from molecu- 
lar clouds interacting with the SNR, that is, CRs with a few 
GeV energies have to reach the molecular clouds. In order for 
Pcut to be smaller than ~1 GeV/c, -Rcsc should be order of 
Ll . Therefore, the effect of the finite-size source considered 
here is important. To estimate the values of those momen- 
tum breaks, we need to specify a model of 7?esc(p). However, 
properties of particle escape depend on models. As we see 
in section 4, when the whole SNR shell interacts with the 
molecular cloud, another momentum break Pbr,i should be 
introduced, which plays a more crucial role than Pcut. 

Given _Reac(p), iesc(p), and Nesc{p), one can obtain 
/ext(i, )", p) and F{p). In section 3 and section 4, we pro- 
vide two models fixing /Jcsc(p), tcsc(p), and Ncsc{p) in the 
context of SNR-molecular cloud interacting system. 



3 MODEL 1: SNR INTERACTING WITH A 
SMALL MOLECULAR CLOUD 

In this section, we consider the case in which a small part of 
a SNR shell interacts with a molecular cloud. In this case, 
the molecular cloud affect neither dynamical evolution of 
the system nor the shock environments to confine CRs, so 
that both Reac (p) and N esc(p) are those for an isolated SNR 
which are obtained from lOhira et al.l (I2OI0I). Here we briefly 
summarize the results of lOhira et al.l (2010|). In the frame- 
work of DSA, CRs are scattered by the turbulent magnetic 
field and go back and forth across the shock front. Once 
CRs reach far upstream where the shock front can not be 
identified as a plane, the CRs can not go back to the shock 
front, that is, the CRs escape from the SNR. In the context 
of DSA, the diffusion length of CR is D{p)/ush, where D{p) 
and Ush are the diffusion coefficient in the vicinity of the 
shock and the shock velocity, respectively. The momentum 
of escaping CRs, Pcsc, is determined by 



:^^^~4sc , (16) 

Msh 

where Iesc is the distance of the escape boundary from the 
shock front. We assume the Bohm like diffusion coefficient 
D — Dop and ^esc = i^Rsh as a geometrical confinement 
condition, wh ere wc adopt k — 0.04 throughout the paper 
(jPtuskin fc Zi rakashvili i2005i ). Then, we obtain 



Pcsc — ^L'q Rsh 



Ush 



(17) 



Although Do has not been understood in detail, pesc starts 
to decrease when the Sedov phase begins at t = tscdov. 



nRshUsh ( Rsh 

Pkn 



/ Rsh Y 

V -Rscdov / 



(19) 



Moreover, we assume that the number of CRs in the momen- 
tum range (nipC, rripC -I- dp) in the SNR is K{Rsh)dp oc R^^^ 
where rup is the proton mass, and that the CR spectrum at 
the shock front is p'". Then, in most of the cases, Ncsc{p) 
has a single power-law form as 



iVcsc(p)0Cp-(^ + ''/"' 



(20) 



The value of /3 dep ends on the inject ion model and ranges 
from -3/4 to 3/13 (jOhira et al.ll2010h . From equation (fTS)). 
we derive 



-Rosc(p) = (1 + K)i?Scdov 



Pkn 



-1/a 



(21) 



Using the Sedov solution and equation (|18[). one finds 

/ \ -5/2c« 

tcsc(p) = iscdov -^— I . (22) 



Pknc 



In this model, Rcsc{p), icsc(p) and Nasc{p) obey single power 
laws, respectively, so that F{p) is simply given by equa- 
tions 112]) or US]). For Pbr,ext > Pbr,2, 



-{1.5S + s + P/a) 



-3/2 



P •■ "'"'Atip)- 

(0.5*+s+(/3-l)/.)^^( )-l/2 

F{p) oc{ ^_(,+^/,) 


For Pbr,2 > Pbr.cxt, 

-(1.5S + s + fi/a)^^,y3/2 



(p > Pbr,Gxt) 
(Pbr.cxt > P > Pbr,2) 
(Pbr,2 > P > Pcut) 
(Pcut > P) 



F(p) OC \ p 





-{a + P/a) 



(P > Pbr,2) 

(Pbr,2 > P > Pcut) 

(Pcut > P) 



(23) 



(24) 



where At(p) = t~ t^scip)- 



4 MODEL 2: SNR EMBEDDED IN A 
MOLECULAR CLOUD 

In this section, we consider the case in which the whole SNR 
shell interacts with the molecular clouds or cavity wall which 
are located at Li < r < L2. Once the CRs encounter the 
molecular clouds, they all are expected to escape from the 
SNR because the high-density neutral gas damps plasma 
waves. In this case, the functional forms of -Resc(p) and 
A^csc(p) are different from those given in the previous sec- 
tion. At an early epoch, the SNR does not interact with 
the molecular clouds and pcsc is again determined by equa- 
tion (|16|l with ^esc = K7?sb . If the shock front gets closer to 
the molecular clouds and Li — R^h < K-R^h, ^csc should be 
replaced with Li — _Rsh. Therefore, 
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Icsc = min [nRsh , Li - R^i,] . 

From equation (|16l) . (|19p and (|25p . we derive 



Pcsc Pkncc 



-Rsh 



^ -RSed. 

The effective power law index of 
dlogpci 



K, \Rsh J 



d log Rsh 



= a + l + 



9csc at i?ah 
1 



K 



(25) 

(26) 
Li/(1 + k) is 

(27) 



-R=h^-£'i/(i+K) 

Since k = 0.04 is small, this index is much larger than a. 
This means that for R^i^ > Li/{1 + k), pcsc rapidly decreases. 
In other words, CRs with p < Phi,i escape at the same time 
at t = tsedov{Li/{l + K)7?sedov)^^^ and Pbr.i is given by 



Pbr.l — PkncG 



Li 



^-RSedov(l +k) 

From equation (|26l) . we find 

(1 + K)i?Scdov (p/Pkncc) 



(28) 



-Rcsc(p) = 



-I/O 



Ll 



(P > Pbr,l) 
(P ^ Pbr.l) 



(29) 



Using the Sedov solution and Eq 

/2q, 
tcsc (p) 



iScdov (p/Pknce) 
iSodov (Pbr.l/Pknoc) 



we find 

(p > Pbr,l) 
(p ^ Pbr.l) 



(30) 



For p < Pbr,i, the spectrum of runaway CRs is the same as 
CR spectrum at the shock because the CRs escape almost 
at the same time. Therefore, in this model, A'esc is approxi- 
mately a broken power law. 



M f ^ ) P-^'^^'"^ (p>Pb.,i) 

iVcsc(p) OC <^ _p/^ -s ( ^ N 

Phr.l P (P^Pbr.l) 



(31) 



Generally speaking, the spectrum is complicated due to 
four breaks. However, in our cases, almost all the CRs can 
reach the molecular cloud, where F{jp) essentially has three 

breaks at Pbr.ext, Pbr,l and Pbr,2. If Pbr,2 < P < Pbr,cxt, 

then F{p)/Nesc{p) has a spectral break due to the break 
of Rd{p)Resc{p) [see equation (I12p ]. Keeping these in mind, 
we derive for Pbr,ext > Pbr.i > Pbr,2 



-(1.5S + s + fi/a), 



-3/2 



P --'At(p)- 

^W«S p-(0.5*+»)Af(p)-V2 



(Pbr.cxt < p) 
(Pbr.l < P < Pbr.cxt) 
(Pbr,2 < P < Pbr,l) 
(P < Pbr,2) 



for Pbr,l > Pbr.cxt > Pbr,2, 

p_(1.55+.+,3/„)^^^p^-3/2 

„. . , ^-(^■^*+^'Af(p)-«''2 



for Pbr,l > Pbr,2 > Pbr.Gxt, 

-(1.5<5+s+,3/a)^^(^p-j-3/2 




-(1.5<5+s) 



Af(p)- 



-3/2 



(Pbr.l < P) 
(pbr.ext < P < Pbr,l) 
(Pbr,2 < P < Pbr.ext) 
(P < Pbr,2) 



(Pbr.l < P) 

(Pbr,2 < P < Pbr.l) 

(P < Pbr,2) 



for pbr,2 > Pbr.cxt > Pbr.l and for pbr,2 > Pbr.l > Pbr.oxt, 

r p-(^-^^+»+''/°'Ai(p)-3/^ (Pbt,2<p) 

Fip)x< p-( = + '9/") (Pbt.l < P < Pbr.2) .(36) 

[ P"" (P<Pbr,l) 

Note that when Pbr,i is small enough, the model 2 is the 
same as model 1 except for pcut (compare equations (23) 
with (35), or (24) with (36)). 



5 GAMMA-RAY SPECTRUM FROM A SNR 

In this section, specifying model parameters, we calcu- 
late the gamma-ray spectrum from a SNR interacting with 
molecular clouds and compare with observed spectra. 

We consider model 2, because observed spectra have 
no low-energy cutoff. From equations Q, (O, (fTi)) . ((28)) . 
H29I) . H30p . and (|3ip . one can calculate the CR momen- 
tum spectrum. For simplicity, we assume that the injection 
model is the therm al leakage model, that is, /3 — 3(3 — s)/2 
(|Ohira et al.ll2010l '). and -Rscdov = 2.1 pc and ^ Sedov = 210 yr 
for al l SNRs. Using the code provided by iKamae et al.l 
(l200a) (see also lKarlsson fc Kamad (120081 ')'). we calculate the 
spectrum of Tr^-decay gamma rays. The normalization of the 
gamma-ray flux is adjusted to flt the data. 

Figures [l] [1 [3] and g] are the results for W51C, "W28, 
W44 and IC 443, respectively. Table 1 is the list of adopted 
parameters and calculated break momentums. We find that 
for all SNRs, the observed gamma-ray breaks can be at- 
tributed to the momentum break at Pbr,i. For W51C and 
IC 443, the model predicts second breaks seen around 
30 GeV and 15 GeV, respectively (see figures [l] and [4]) . For 
W51C, the break comes from the momentum break at Pbr,2 
that is the effect of a finite-size emission zone, and for IC 443, 
the break comes from the momentum break at pbr.ext that 
is the effect of a finite-size source. Note t hat Pbr,2 was also 
predicted bv lAharonian fc AtovanI (| 19961 ). while other two 
breaks , Pbr.i and Pbr.ext, are first introduced in this paper. 
For W28 and W44 (figures [5] and [S]) , Pbr,2 and Pbr.oxt are 
much smaller than the threshold of the tt*^ production, so 
that the spectrum is the same as that of the point source 
solution (equation [2}. Our model seems consistent with ob- 
vsfetvations for all the SNRs. The observed data are fitted 
with the same parameters concerning the escape and the 
acceleration, a, /3, k, and s. Hence, all the SNRs could have 
the same source spectrum above pbr,i(~ 1 — 10 GeV), re- 
quired for explaining Galactic CRs. The observed diversity 
of the gamma-ray spectra comes from five parameters re- 
(33) lated to CR propagation around SNRs (x and 5) and/or 
environments (Li, L2, and tagc). 



(34) 



for Pbr.oxt > Pbr,2 > Pbr 

-(1.5S+a+P/a). 



-3/2 



F(p) OC ^ ^_(,+^/,) 



p-'-—^'^^'-'At{p)- 

(0.5S+s+(P-l)/a)^^l^pyl/2 



(Pbr.cxt < P) 
(Pbr,2 < P < Pbr.cxt) 
(Pbr.l < P < Pbr, 2) 
(P < Pbr,l) 



6 DISCUSSION 

For the four SNRs detected by Fermi, we have shown that 
they can be explained by an accelerated-CR spectrum with 
s ~ 2.2 and the corresponding runaway-CR spectrum with 
s + PI a ~ 2.38. It is found that s = 2.1 — 2.4 also gives ac- 
ceptable fits. Those spectral indices are consistent with those 
/^,the source spe ctrum of Galactic CRs, wh i ch is implied 
)in GALPROP JStrong fc Moskalenkd[l99i: IStrong et al " 



2OOOI ) and others (e.g.. IShibata et al.l 12009 : iPutze e~ 
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Table 1. Model parameters and characteristic momentums 



SNR 



(pc) 



L2 
(pc) 



f-agc 

(kyr) 



Pbr.l 

(GeV/c) 



Phi, 2 
(GeV/c) 



Pbr,cxt 

(GeV/c) 



W51C 


6.5 


1.2 


0.04 


2.2 


0.22 


0.1 


14.7 


23.7 


31.5 


13.1 


124 


21.6 


W28 


6.5 


1.2 


0.04 


2.2 


0.19 


0.9 


18.5 


26.9 


63.0 


2.96 


< 1.43 


< 1.43 


W44 


6.5 


1.2 


0.04 


2.2 


0.40 


1.0 


12.4 


16.2 


23.1 


39.8 


< 1.43 


< 1.43 


IC443 


6.5 


1.2 


0.04 


2.2 


0.62 


0.01 


11.5 


14.7 


23.1 


62.9 


2.93 


151 




10' 10^ 

Energy [MeV] 




10"* 10= 10° 

Energy [MeV] 



Figur e 1. Compar ison of the model results (s olid line) with Ferm i 
(red) l lAbdo et al.| [2009b) and HESS (blue) l lFiasson et al.ll2009t) 
observations for SNR W51C. 




10' 10° 

Energy [MeV] 



Figur e 2. Compar i son of t he model results (soli d line) with Fermi 
(red) llAbdo et al.l l2010bl) and HESS (blue) llAharonian et al.l 
|200^ observations for the source N of SNR W28. 



l2009r ). within uncertainties in the CR propagation. The in- 
ferred values are somewhat steeper than that predicted by 
the classical DSA (s = 2), and several mecha nisms have 
been proposed to explain such steeper spectradKirk et al.l 
19961: IZirakashvili fc PtuskinI l2008l : lOhira et al.l l2009bl ld: 



Fuiita et a l. 2009). However, our results are based on several 



approximations, and we have investigated only four middle- 
aged SNRs, so that results on indices are not conclusive yet. 
Further theoretical and observational studies are obviously 
required in order to discuss underlying acceleration mecha- 



Figure 3. Comparison of the model results (s olid line) with 
Ferm i (red) llAbdo et al.ll2010d), Whipple (gre en) llSucklev et al.l 
Il998l). HEGRA (blu e) llAharonian etal.|[20o3) . Milagro (purple) 
JAbdo et ai] l2009a^ observations for SNR W44. The data of 
Whipple, HEGRA and Milagro are upper limits. 
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Figur e 4. Comparison of the model results (solid lin e) with Fermi 
(red) llAbdo et al.ll2010a^ . MAGIC (e reen squares) JAlbert et al.l 
l2007t ) and VERITAS (purple circles) l|Acciari et al.ll2009h obser- 
vations for SNR IC 443. 



nisms and reveal the link between those middle-aged SNRs 
and observed Galactic CRs. 

In order to fit the data, the parameters concerning the 
diffusion coefficient of CRs just around the SNR, x ^^'^ "^i 
are different from galactic mean values, x ~ 1 ''■nd ^ ~ 0-5 
(see Table 1). Possible explanation of these discrepancies is 
that the runaway CRs amplify the magnetic fluctuation, so 
that the diffusion coefficient in the medium with the high 
CR density can be different from the galactic average. The 



6 Y. Ohira et al. 



higher the CR density is, the more turbulent magnetic field, 
and the more fi'equent the scattering of CRs, resulting the 
smaller dilTusion coefiicient. On the other hand, the mag- 
netic turbulence becomes weak if the neutral particles im- 
pedes the magnetic-field amplification. Further studies are 
needed to understand physics of partially ionized interstellar 
matter illuminated by CRs. 

Some important processes neglected so far should 
be mentioned. At first, in § [J] the diffusion coeffi- 
cient of CRs in the interstellar medium, -Dism(p), is 
assumed to be sp atially uniform and time-independent. 
iFuiita et al.l (|20ld ) demonstrated that CRs that have es- 
caped from an SN R excite plasma waves and reduce the 
diffusion coefficient jWentzejl 19691 : iKulsrud fc Pearcelll969l : 
IPtuskin et al.ll200g ). which is responsible for non-uniform, 
time-dependent diffusion coefficient. Secondly, in § [3] and 
§ m we assume the diffusion coefiicient just around the 
shock has a form of Eq. (|19|) . Although the magnetic 
field amplifica tion aroun d the shock is studied by the lin- 
ear an alysis l|Belll l2004l : iReville et all l2007l: lOhira et all 
,2009bl) and many numerical simulations (Lucek & Bell 200(3: 
Giacalonc & Jokipii 2007; Nicmicc ct al. 2008; Rcvillc ct al.l 



20081: iRiauelme fc Spitkovskvl |2009|: llrioue et all 12009. . 
Ohira et al.ll2009al : IVladimirov et ani2009l 'l. the saturation 



of the magnetic-field amplification and the diffusion coeffi- 
cient has not yet been understood in detail. Finally, in § [3] 
and § m we assume that the CR spectrum at the shock has 
a single power law, p^" , at any time. However, the situation 
may not be so simple. The shock compression ratio decreases 
as the shock velocity decre ases. As a result, the spectrum 
becomes steeper with time (jFuiita et al.ll2009l ). In addition, 
neutral particles have great i mpacts on the mag netic field 
and the particle acceleration (jOhira et al.ll2009bl la). More- 
over, nonlinear DSA should be c onsidered when CRs affect 
the s hock structure JDrurv fc Vo lk 1981: iMalkov fc Drurvl 
[2OOIJ), along with properti es of magnetic turbulen ce in mag- 
netic field amplification (|Vladimirov et al.l |2009| ). Further 
studies are needed to clarify the significance of these effects. 
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